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Abstract. Among the key features of hot and dense QCD matter produced in ultra- 
rclativistic heavy-ion collisions at RHIC is its very low shear viscosity, indicative of 
the properties of a near-ideal fluid, and a large opacity demonstrated by jet energy 
loss measurements. In this work, we utilize a microscopic transport model based on 
the Boltzmann equation with quark and gluon degrees of freedom and cross sections 
calculated from perturbative Quantum Chromodynamics to simulate an ideal Quark- 
Gluon-Plasma in full thermal and chemical equilibrium. We then use the Kubo 
formalism to calculate the shear viscosity to entropy density ratio of the medium 
as a function of temperature and system composition. One of our key results is that 
the shear viscosity over entropy-density ratio rj/s becomes invariant to the chemical 
composition of the system when plotted as a function of energy-density instead of 
temperature. 
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1. Introduction 

Ultrarelativistic heavy ion collisions at the Relativistic Heavy Ion Collider (RHIC) are 
thought to have produced a Quark Gluon Plasma (QGP) with the characteristics of 
a near ideal fluid|U El [3j H]. One of the most important current challenges in QGP 
research is to quantify the transport coefficients of this novel state of matter. Recently, 
attention in the field has been primarily focused on the shear viscosity to entropy density 
ratio r]/s. Calculations utilizing certain strongly coupled supersymmetric gauge theories 
with gravity duals j5] postulate a lower bound of 77 min = s/4tc for this quantity, often 
referred to as the KSS bound [6]. Relativistic viscous hydrodynamic calculations have 
confirmed very low values of rj/s in order to reproduce the RHIC elliptic flow (t^) data 
[H[8j|9]. However, current calculations assume a fixed value of rj/s throughout the entire 
evolution of the system and neglect its temperature dependence, which anyhow cannot 
be obtained from the AdS/QCD calculations. 

It has been argued on very general grounds that r]/s should exhibit a minimum 
at the phase-transition and should rise for decreasing temperature in the hadronic 
phase and for increasing temperature in the deconfined phase [ID] . Recent calculations 
of rj/s using microscopic transport theory in the hadronic sector have confirmed this 
expectation and have provided quantitative guidance with respect to the temperature 
dependence of rj/s in and out of chemical equilibrium [TT]. At temperatures above Tc 
there remains a large uncertainty regarding the value and temperature-dependence of 
rj/s: in the limit of very high temperatures, the interactions between the constituents 
of QCD matter should be calculable with perturbative methods [12] , whereas in the 
strongly interacting regime near Tq rj/s could be close to the KSS bound. It is our goal 
in this paper to explore the temperature dependence of QCD matter for temperatures 
above Tc at which quasi-particle degrees of freedom become viable. For our calculation 
we will use a microscopic transport model, the BMS implementation of the Parton 
Cascade Model (PCM) [13], which has the advantage that it can describe both, out 
of equilibrium QCD matter created in collisions of ultra-relativistic heavy-ions, as well 
as equilibrated QCD matter at a fixed temperature. For the calculation of the shear 
viscosity we shall use the Kubo formalism utilized in our previous calculation of r] and 
rj/s in the hadronic sector [TT] . 

2. The Parton Cascade Model 

The Parton Cascade Model (PCM) [TTJ [15] is a microscopic transport model which 
is used to simulate the time evolution of a system of quarks and gluons utilizing 
the Boltzmann equation. Here, we use the BMS implementation, which has been 
applied successfully to the calculation of direct photon production [161 E] an d baryon 
stopping [18] at RHIC. PCM implementations by other groups have been used to 
address collective flow, parton energy-loss as well as multi-particle effects on transport 
coefficients [19l [2Ql ETJ [22l [23l EH [25l [26]. Our calculation focuses on an ideal Quark- 
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Gluon-Plasma, i.e. a gas of u, d and s quarks and anti-quarks as well as gluons at 
fixed temperature T in full thermal and chemical equilibrium. In addition, we also 
study a one-component gluon plasma in thermal and chemical equilibrium. For our 
transport calculation we define a box with periodic boundary conditions (to simulate 
infinite matter) and sample thermal quark and gluon distribution functions to generate 
an ensemble of particles at a given temperature and zero chemical potential. 

The time evolution of our system is described by a Boltzmann transport equation: 

„ d 



P 



dx^ 



F k (x,p) 



(i) 



processes: i 

where the collision terms are nonlinear functionals of the phase-space distribution 
function. Although the collision term, in principle, includes factors encoding the Bose- 
Einstein or Fermi-Dirac statistics of the partons, we neglect those effects here. 
The collision integrals have the form: 

(27T) 4 
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Si is a statistical factor defined as ^ = T7 K™\ K^l with fQ n,out identical partons of 

species a in the initial or final state of the process i, excluding the kth parton. 

The matrix elements for the full Quark-Gluon-Plasma calculation \M\ 2 account for 
the following processes (note that all t and u-channel processes are included as well): 







The corresponding scattering cross sections are expressed in terms of spin- and 
color-averaged amplitudes \M.\ 2 : 

In the case of a pure Gluon Plasma (GP), we use the following differential cross 
section in order to facilitate comparisons with PCM implementations by other groups 



d a 99^99 2 9 1 
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For the transport calculation we also need the total cross section as a function of s 
which can be obtained from (151): 



MS) = £ / 773 dQ 2 . (7) 




The amplitudes for the above processes have been calculated in refs. [27] [28] for 
massless quarks. Note, however, that for the kinematics of the system, the light quark 
masses are treated explicitly. We regularize the IR divergence of the cross sections with 
a temperature dependent Debye mass mr>. We shall use two different expressions for 
vtld - the first one is a Debye- mass for particles obeying Boltzmann statistics which has 
been used in [211 E3J EH [251 EE] and which we utilize to allow our results to be compared 
to these calculations: 




m D {T) = \—a s T* , (8) 



and the second one is the standard Debye mass used in pQCD calculations for systems 
of quarks and gluons: 



m D (T)=gT^/(2N c + N f )/6 (9) 

The first mj) parametrization we shall refer to as Boltzmann-m^ whereas the second 
parametrization we shall refer to as regular mjj. In both cases the coupling constant is 
defined as 

a ° = £ (10) 
and can either be chosen as a constant parameter, or to have the following temperature 
dependence [10J: 

g 2 8tt 2 \A t J 9vr 2 V \A T j) 
with At = 30 MeV. 

3. Methodology 

3.1. Shear-viscosity 

In order to calculate the shear-viscosity of our system, we employ the Kubo-formalism 
[29| [30] . The methodology for applying the Kubo-formalism to infinite QCD matter 
modeled by microscopic transport theory was developed in [31J and has already been 
successfully applied to calculating the shear viscosity for a hadron gas [311 [11]. The 
momenta of all the partons in our system are tracked over the course of 20 time-steps 
of 0.5 fm/c. Knowing the momenta p of all partons in the system, the discretized 
stress-energy tensor of our system can be calculated at each timestep: 

71 = ldp—^f(x,p)^ -XJ ~ 
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Figure 1. Kubo correlators for a system of gluons at various temperatures. The 
correlators show a clear exponential decay. 



The stress-energy tensor now allows for the construction of the stress energy tensor 
correlator used in the Kubo formalism [291 [30J to calculate the shear viscosity via 

1 f poo 

r ] = -jd 3 rj o dt{n*»(0, 0)^(f, t)) eqml (13) 

where the average is taken over many events and T is the temperature of the system. 
A selection of correlators for different temperatures are shown in figure 1. Based on the 
observed exponential decay of the correlators, we fit these with the following analytic 
function: 

(ir xy (0,0)ir xy {r,t)) ~ (14) 

with the parameter tr called the relaxation time of the system. Inserting this function 
into the expression for r\ yields: 

17= ^7^(0,0)7^(0,0))^ (15) 

with the volume V and temperature T being input parameters of our calculation and 
tr and {n xy (0, 0)ir xy (0, 0)) determined by our exponential fit (see figure 1). 

3.2. Entropy 

The entropy-density of our system is being calculated via the Gibbs entropy density: 

s = y (16) 

Since the system is initialaized with net baryon density, p = 0, the pressure P and 
energy- density e are the only quantities which need to be extracted from our simulation 
according to 




Figure 2. Left: entropy as a function of temperature for a Gluon-Plasma and a three 
flavor Quark-Gluon-Plasma. The symbols denote the entropy of the system evaluated 
using the Gibbs formula and the lines represent the Stefan Boltzmann expression for 
the entropy density with v(T) — 16, and Quark-Gluon-Plasma, v(T) = 47.5 Right: 
shear-viscosity 77 as a function of temperature for the Gluon Plasma (GP) and the 
Quark-Gluon-Plasma (QGP). 



The left frame of figure 2 shows the resulting entropy-density for a Gluon Plasma (GP) 
as well as a Quark-Gluon-Plasma (QGP) containing gluons and three light quark flavors. 

We can compare the entropy density directly obtained from our calculation to the 
Stefan Boltzmann entropy density given by: 

a = ^(T)T 3 (18) 

where u(T) = N^+^Nf and Nf, and Nf are the bosonic and fermionic degrees of freedom. 
The results of the comparison are shown in the left frame of and figure 2 and show good 
agreement between the Stefan Boltzmann entropy and our system entropy calculated 
via the Gibbs relation. 



4. Results and Discussion 



We will present our results predominantly for two modes of calculation: the first 
one, denoted by GP, is for a gluon plasma with a regular Deybe screening mass and 
temperature dependent coupling. The second mode, denoted by QGP, uses the same 
temperature-dependent parametrizations for mo and the coupling constant, but is for 
a quark-gluon plasma with three light quark flavors. We consider the QGP mode to be 
the most realistic mode of calculation presented here. 

The right frame of Figure 2 displays the results of our calculation for the shear 
viscosity for a system of gluons as well as a system of gluons and three light quark 
flavors. Both, the shear viscosity, as well as the entropy density depicted in the left 
frame of figure 2, rise strongly as a function of increasing temperature. It is interesting 
to note that the QGP shows a significantly higher shear- viscosity than the GP for a given 
temperature, which seems counter-intuitive given the larger particle density of the QGP, 
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Figure 3. Left: rj/s as a function of temperature for a Gluon-Plasma (full circles) 
and Quark-Gluon-Plasma (full triangles), compared to a HTL calculation (solid line) 
[12]. Right: same result, but plotted vs. energy density. 



but is probably due to the smaller interaction cross-sections among the quarks of the 
system. 

Having calculated both, the shear-viscosity as well as the entropy-density of our 
system, we can now turn to the ratio rj/s, made famous by the KSS bound: the left frame 
of figure 3 shows 77/s as a function of temperature for the GP (full circles) and QGP 
(full triangles) case, compared to an analytic HTL calculation of a three-flavor QGP 
[12] (solid line). The QGP calculation of 77/s shows a monotonous rise as a function of 
temperature with a slope very similar to that of the HTL calculation. The differences in 
absolute value between the two can be understood by considering the NLL corrections 
present in the HTL calculation into account, which are missing in our simulation of QCD 
matter. The Gluon Plasma exhibits an upturn in 77/s for temperatures below 500 MeV 
- we attribute this unexpected rise towards lower temperatures to a breakdown in the 
perturbative approximations present in our calculations, emphasized by taking the ratio 
of two quantities, which both, when calculated and viewed separately (see figure 2), seem 
to exhibit a smooth behavior in the low temperature domain below 500 MeV. 

Comparing 77/s of a GP and a QGP at the same temperature may be misleading due 
to the significantly larger parton density present in a QGP. Therefore in the right frame 
of figure 3 we compare rj/s for the two scenarios at equivalent energy- density and find for 
energy-densities above 35 GeV/fm 3 excellent agreement between the two systems. This 
is of particular relevance since the flavor composition of the deconfined QCD matter 
created in ultra-relativistic heavy-ion collisions is by no means fully established and 
may change strongly as a function of time - from a gluon-dominated system being 
created by the decay of a Color-Glass-Condensate to a QGP in full thermal and chemical 
equilibrium as the system progresses in its hydrodynamic expansion. Our result indicates 
that rj/s, a quantity which controls the hydrodynamic evolution of the system, should 
be fairly robust with respect to its flavor composition when taken as a function of 
energy-density instead of temperature. 
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Figure 4. Left: effect of different Debye mass parametrizations and temperature- 
dependent coupling on rj/s. Right: Effect of coupling constant on rj/s. 



In the left frame of figure 4 we study the effects of our different Debye-mass 
parametrizations and the temperature-dependence on rj/s. For this purpose we calcuate 
rj/s for a gas of gluons with a Boltzmann Debye-mass at fixed coupling of a s = 0.3 (i.e. 
the same system as e.g. in J2TJ [22j [23j |2U [25] ) and compare this calculation to one 
with our default parametrization (taken with N c = 3 and Nf = 0). We find the effect 
of the different mp parametrizations to be small, on the 10% to 15% level, with the 
Boltzmann Debye-mass giving systematically smaller values of rj/s. If we now replace 
the fixed coupling with a temperature dependent coupling constant, the value of rj/s 
increases roughly by a factor of 2. 

So far we have restricted our investigation to a purely perturbative partonic 
system with the respective temperature-dependent screening masses and coupling 
strengths. However, there are strong indications that the medium created in ultra- 
relativistic heavy-ion collisions is non-perturbative in nature - at the very least within 
the temperature range from Tq to approximately 3 — 4Tc which is covered by our 
calculations. One method to explore the behavior of rj/s at stronger coupling is to treat 
the coupling constant in our calculations as a free parameter and then study rj/s at fixed 
temperature as a function of the coupling constant. The right frame of figure 4 shows 
rj/s as a function of coupling strength for the gluon plasma and the quark-gluon-plasma. 
In the strong coupling limit, in particular for the gluon plasma, values of rj/s <C 1 can 
be obtained, yielding results compatible with a fluid-dynamical analysis of RHIC data 
[9]. Similar results have been obtained by [20], who directly increased the gluon-gluon 
scattering cross-section, which is equivalent to an increase in the coupling constant 
or by directly studying rj/s as a function of the coupling constant [22]. One should 
note, however, that treating the coupling constant as a free parameter yields a medium 
with inconsistent characteristics, since the particle density and screening mass are still 
controlled by the temperature, whereas the coupling is not. Also, for large couplings 
a s ~ 1 the perturbative assumptions underlying the PCM are not valid anymore. 

The final question we wish to address is the effect the angular distribution of 
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Figure 5. Coupling constant and angular distribution dependence of rj/s for a gluon 
plasma: The choice of a fixed coupling constant affects r\/s to a far larger degree than 
the transition from the regular angular distribution to isotropic scattering. 



the scattering partons , i.e. the exact form of the differential matrix-element, has 
on rj/s: figure 5 shows a comparison between our default gluon plasma calculation 
compared with results at fixed coupling and regular matrix elements as well as with 
the same matrix elements using an isotropic angular distribution. The main effect, a 
reduction of rj/s by a factor of 3, comes from the choice of a fixed coupling constant 
vs. a temperature dependent one. Changing from forward-backward peaked to isotropic 
scattering provides an additional 10% - 20% effect, but not a dramatic reduction in rj/s. 
This finding is of particular interest in the context of work done by [2U [221 [231 GU E5] , 
employing a PCM including radiative corrections (i.e. 2—7-3 and 3—7-2 processes) 
and an implementation of the LPM effect, which manifests itself in a near- isotropic 
angular distribution for the third particle in the outgoing channel. It has been shown 
that the isotropic emission of the radiated particle in this implementation of the LPM 
effect plays an important role for the thermalization of the system. Our results utilizing 
solely 2—7-2 scattering processes indicate that it is most likely a combination of multiple 
effects - choice of fixed coupling constant, the isotropic angular distribution of the LPM 
implementation, and the inclusion of radiative corrections, which is responsible for the 
low viscosity findings of these calculations [22| [23] . 

5. Summary 

We have utilized the parton cascade model to simulate a perturbative quark-gluon- 
plasma in full thermal equilibrium and have extracted its shear- viscosity as well as the 
shear-viscosity over entropy-density ratio as a function of temperature. We find that 
our results depend significantly on the details of the calculation, i.e. choice of coupling 
constant and parametrization of the Debye screening mass as well as the degrees of 
freedom (gluon vs. quark-gluon plasma). One of our key results is that the shear 
viscosity over entropy-density ratio rj/s becomes invariant to the chemical composition 
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of the system when plotted as a function of energy- density instead of temperature. The 
values we obtain for rj/s are higher than those expected for the near ideal fluid observed 
at RHIC; in particular they are not compatible with rj/s values extracted from viscous 
fluid dynamics analysis [9] of elliptic flow data. By increasing the coupling constant 
we find values of rj/s that are compatible with the RHIC data, but only for values of 
the coupling at which the perturbative assumptions of the PCM may not anymore be 
valid. Inclusion of quantum- coherence effects, such as the LPM effect, multi-particle 
scattering processes or turbulent color fields leading to an anomalous viscosity [32l 133] 
may explain the origin of the observed small 77/s values, but the final determination of 
this question remains to be settled in future work. 
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